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■ Summary. The origin of the abundance discrepancy in nebulae is one of the key 
problems in the physics of photoionized nebulae. In this work we have analized a 
sample of Galactic and extragalactic H II regions where the abundance discrepancy 
have been measured, and we discuss the two main scenarios proposed to explain 
such discrepancy: temperature fluctuations over the observed volume of the nebulae 

^ . and chemically inhomogeneous inclusions. 
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^ [ 1 Overview 

H, 

\ Almost forty years ago, [15] claimed that in the presence of temperature vari- 

■ ations over the observed volume the abundances derived for nebulae from 
the analysis of intensity ratios of colisionally excited lines (CELs) were un- 

O ^ derestimated if these temperature fluctuations were not considered. This is 

' | due to the strong dependence of the intensity of those lines on the assumed 

; I ■ electron temperature. This possibility is puzzling, considering that the anal- 

\ ysis of CELs is the standard method for deriving ionic abundances in ionized 

■ nebulae. On the other hand, intensity ratios of recombination lines (RLs) are 
almost independent of the temperature structure of the nebula, and could, in 
principle, be more appropiate to derive the "real" abundances of the nebulae. 

rS | Measurements of RLs on H II regions and planetary nebulae (PNe) in the last 

■ years have found that abundance determinations from RLs are sistematically 
larger than those obtained using CELs, independently of the ion considered 
(see e.g. [10], and [6] and references therein). 

[20] have presented an scenario to explain the abundance discrepancy in 
Hii regions, which is based on the implicit assumption on the existence of 
temperature fluctuations produced by chemical inhomogcncitics in the ionized 
gas. These authors proposed a model for 30 Doradus, in which they postulate 
the existence of a low temperature, metal-rich gas embedded in an "ambient" 
medium with lower density and metallicity and higher temperature. According 
to [20], this denser component comes from Type II SN ejecta, that has not 
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been mixed with the interstellar medium and is in pressure equilibrium with 
the "normal" chemical composition gas. This droplets would be responsible of 
most of the RLs emission, and due to their low temperature they would not 
emit CELs. 

In this paper we present results on deep VLT high resolution spectropho- 
tometriy obtained by our group in the last years ([1, 8, 5, 6, 7, 13]) to discuss 
the validity of the two scenarios described above to explain the abundance 
discrepancy problem. 

2 Correlations of the ADF with nebular properties 

The different scenarios proposed to explain the abundance discrepancy predict 
a different behaviour of the ADF with respect to different nebular properties. 

The metallicity dependence of the t 2 parameter was suggested by [3] , who 
found that photoionization models could reproduce temperature fluctuations 
similar to those observed in the case that electron temperature is less than 
~9000 K, raising on colder -more metal-rich- nebulae. This result does not 
agree with the observations in both H II regions and PNe. In figure 1 we show 
ADF vs. ++ /H + and excitation, for all the Galactic and extragalactic Hii 
regions for which ++ abundances from O II RLs are available. Despite the 
narrow range of parameters covered, we have found no correlation between 
ADF(0 ++ ) and the rest of parameters. 




Fig. 1. ADF(0 ++ ) vs. excitation (left) and vs. 0++/H+ ratio (right). Filled circles: 
Galactic Hn regions; rest of symbols: extragalactic Hn regions from the literature 
(see [13] and references therein). 

On the other hand, several authors have found that the large range of 
ADFs found for PN is because the ADF increases monotonically with metal- 
licity (see e.g. [12]). In the case of Hn regions, this effect has not been found 
(see Figure 2), at least in the studied metallicity range. We have also rep- 
resented ADF vs. T e (high) /T e (low) and we have found no correlation (see 
Figure 2), indicating that the large scale variations of T e due to the natural 
gradients of this parameter along the nebulae seem not to be related with 
the ADF. Finally, we have also check that the ADF do not depend on the 
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assumed temperature (see Figure 2), which discards systematical effects in 
the abundance determination from CELs. 
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Fig. 2. ADF(0 ++ ) vs. oxygen total abundances from CELs and RLs (left) and vs. 
T e (high) and T e (low)/T e (high). Symbols are the same than in Figure 1. 

It is clear that it is necessary to increase the sample of H n regions in which 
RLs have been measured with a good signal-to-noise. A systematic search of 
RLs in extragalctic H n regions would permit us to cover a larger fraction of 
the nebular volume, to increase the metallicity range. 

Finally, [17] proposed that photoelectric heating due to dust grains could 
increase T e in the zones near the central star of a PN. [16] found that the 
ADF was not correlated with the effective temperature of the ionizing star 
for a sample of PN. For our sample of Hn regions we have not found a cor- 
relation between the ADF and the effective temperature of the ionizing star 
-or spectral type, which could give an idea of the radiation hardness- of each 
nebulae (see [4]). 




3 Implications for the temperature fluctuations scenario. 

Many of the arguments against temperature fluctuations are based on that 
photoionization models can not reproduce observed i 2 's, but photoionization 
models could not be enough realistic, and additional energy sources could be 
necessary to explain the observed discrepancies among observational data and 
photoionization models predictions (see [21] and references therein). 

In the last years, several observational evidences have been presented 
against temperature fluctuations in PN (see [10] and references therein). On 
the other hand, very recently, [13] reported some reasons against the two- 
phase model of [20]. Here, we present additional objections to the two-phase 
model. By definition, under the temperature fluctuations scenario, the ADF 
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should be related to the excitation energy (see [15]); on the other hand, un- 
der the density fluctuations scenario ([22]) the abundance discrepancy should 
be higher if the abundance has been computed from a CEL with an upper 
level with low critical density, n EjCr , so ADF and critical density should be 
inversely proportional. [12] and [11] compared several abundance determina- 
tions in PNe (CELs in UV, optical and far-IR), and showed that the ADF 
was not related with the excitation energy, nor with n e ^ cr . As it has been 
pointed out in previous works (see [6, 7]), we can not compare abundances 
derived from CELs in different ranges (UV, optical or far-IR) of the same ion, 
because UV and far-IR observations covering the same slit areas in the optical 
range are not available. Nevertheless, we can compare the ADFs obtained for 
Galactic and extragalactic Hii regions, with the excitation energy and with 
the critical density, n e ^ cr , of the upper level of the main CEL of each ion. 
In Figure 3 we shown a possible correlation (r=0.6) between the ADF and 
the excitation energy, E ea; ; there is also an apparent correlation (r=0.5) with 
n e;Cr , contrary to the predictions of the density fluctuations theory. These 
results are an additional probe that the abundance discrepancy scenario in 
H II regions should be different to those proposed for PNe, so results on PNe 
can not be extrapolated to Hn regions. 
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Fig. 3. ADF for different ions vs. critical density (left panel) and vs. excitation 
energy of the upper level of the line (right panel) . 



3.1 Temperatures from On RLs. 

The chemically inhomogeneous or two phase model ([11], [18] and [20]) pre- 
dicts T e (RLs) < T e (Hei) < T e (H i) < T e (CELs) ([9]), with the difference 
between temperatures being proportional to ADF (see e.g. Figure 8 of [11]). 

[23] used -for the first time- the temperature dependent ratio /(A4089.29) / 
/(A4649.14) to derive T e of the ionized gas from which O n lines arise in PNe 
(T e (RL)). These authors found very low T e in two H-deficient knots of the PN 
Abell 30. Later, in subsequent papers they have found similar results in other 
PNe (e.g. [19, 24]). This method present several difficulties that were pointed 
out by [19]. 
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The determination of T e from O n RLs in H n regions is very difficult 
because the lines are usually much fainter than in bright PNe. fn Figure 4 we 
have represented the dependence of the I(A4089.29)//(A4649.14) ratio with 
T e , and we have superimposed the ratios obtained in NGC 3576 and the Orion 
Nebula, the only objects in our sample in which it has been possible to detect 
and measure the O n A4089.29 line. Moreover, we have represented the line 
ratios measured in the giant Hn region 30 Doradus [14] and in several PNe 
from the literature: (most of them from the sample of [19] and a sub-sample 
with ADF >4 from [24]). 




Fig. 4. On A4089.29/A4649.14 ratio as a function of electron temperature. Solid 
line in the three panels represents the theoretical ratio. Observations are shown as 
error boxes. Vertical arrow shows an object with a ratio outside the scale. Panels 
are described in text. 



As it is shown in the upper-left panel of Figure 4, PNe with large ADFs 
show, in general, very low T e (On), which is consistent with the two-phase 
model. On the other hand, PN with moderate values of the ADF (upper-right 
panel of Figure 4) show T e (On) that are, in most of the cases, smaller than 
T e (CELs) and T e ([Hi]). In the case of Galactic Hn regions, the temperature 
obtained for the Orion Nebula is somewhat larger than that obtained from 
CELs; the ratios obtained for NGC 3576 and 30 Doradus are higher, but the 
Oil A4089.29 line in both nebulae is affected by charge transfer effects in the 
CCD (due to the presence of very bright emission lines in adjacent orders) that 
can not be subtracted due to their low signal-to-noise ratio (see [4]), With the 
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aim of obtain additional data for H II regions, we have represented the ratios 
for two additional positions in the Orion nebula: the position labelled by [2] 
as "Orion 1", and a slit position placed on the Orion bar (24" N and 12" O 
from # 2 Ori A, labelled as "Orion Bar"). It can be see that T e values for this 
regions are similar to those obtained from the values obtained by [1]. 

4 Conclusions. 

We have presented here the results from deep high resolution spectropho- 
tometry of H II regions. From the analysis of these data, it is clear that the 
abundance discrepancies in H II regions do not show the same behavior than 
in PNe, and could be explained assuming temperature fluctuations over the 
observed volume of the nebula. Moreover, the results obtained for PNe can 
not be extrapolated for Hil regions; in fact, it seems that the two-phase model 
scenario proposed by [20] is not suitable to explain the abundance discrepancy 
in Hn regions. 
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